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ABSTRACT 

The fundamental mode of tadial pulsation has been calculated 
for convective envelopes in adiabatic equilibrium including 
radiation pressure. The theoretical Q values are compared with 
the observational Q values of variable M supergiants in galactic 
associations. The variability of these massive red supergiants 
can be explained on the pulsation hypothesis if a suitable 
readjustment of some of the observed periods is made. Pulsation 
seems to be suggested by (1) the cyclical variations in magnitude, 
spectrum, and radius in individual stars, and (2) a period- 
luminosity relation and (3) a period-spectrum relation for 


the group of stars. 
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I . INTRODUCTION 

It has long been assumed that the variability of bright red 
supergiants associated with young populations is due to radial 
pulsation. Some photometric and spectroscopic studies have 
strongly supported the hypothesis. In the interest of bringing 
further evidence to bear on the point, this paper presents a 
theoretical and observational determination of the period- root-raean- 
density quantity 

— L 

Q = Per (p/p^)^ 

for masbive stars. The theoretical Section II contains a brief, 
but complete, discussion of the fundamental pulsational properties 
of convective envelopes in adiabatic equilibrium, while Section 
III compares these results with the observations pf variable M 
supergiants in galactic associations, whose distances are known. 

II. THEORETICi'L MODEL ENVELOPES 

The structure and pulsational characteristics of convective 
stellar envelopes are very difficult to calculate in detail. , 
However, the superadiabatic region near the surface is, in general, 
rather small and the solutions converge rapidly toward an adiabatic 
solution in the interior. The particular adiabat which is attained 
depends on the properties of the superadiabatic region and may be 
characterized by a parameter K (Hayashi, Hoshi, and Sugimoto 1962) . 
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in the following, we make two simplifying approximations. 

First, we assume adiabaticity throughout the envelope and leave 
K a parameter to be specified. Second, we consider the equation 
of state to be represented simply by a mixture of perfect gas and 
radiation. The latter approximation is certainly valid in the 
interior. Previous work on the pulsational characteristics of 
completely adiabatic envelopes neglected radiation pressure. It 
is necessary to include it here explicitly because of its importance 
at the high masses which are typical of red supergiants. Since 
the few equations required here are scattered in five different 
publications, it is useful and convenient to assemble and discuss 
them briefly together. Defining l-P as the ratio of radiation 
pressure to total pressure, we may describe the state of the gas 
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by the two adiabatic exponents (Chandrasekhar 1939) 

^ 32-24f?-3P^ _ 32-243-33^ ,,, 

”l 24-21^ ' ^2 24-18 3-3 32’* 


The basic differential equation^ describing the structure of 
a convective envelope in adiabatic equilibrium may be expressed in 
non-dimensional form, 


q dx qt * 


X dp ^ 3q - 
p dx xt 


V, 



p dt “T 2 -I 


= n+1, 


( 2 ) 



5 t 


where we have adopted the usual Sdhwarzschild variables and have 
included radiation pressure (e.g., Hayashi ^ 1962). We shall 

assume zero boundary conditions at the surface, p = t ® 0 at 
q = X = 1. Substituting the expression for the second adiabatic 
exponent, we integrate the last equation to obtain 

p = Et5/2(i + (3) 

where E is the constant of integration* Near the surface we have 
approximately 

q = 1/ t = .J^(^l) , V = (4) 

A convective envelope is specified by the parameters E and 
A, where 

E = 4TTG'’^^(^H/k) A = ( 4n/3) (pH/k) V . (5) 

A solution loops in the U-V plane toward the point U = 0, 

V = (n+l)g, where (n+l)^ = 2.5 if A = 0 and (n+l)g = 4 otherwise. 
In the latter case, the relative radiation pressure 1-3 increases 
monotonically. from zero at the surface to unity at the asymptotic 
limit of mass fraction, q . For this reason, we define a mean 
value, = J(l-3)dx, which indicates the relative importance of 

radiation pressure. As E increases with A held fixed, the mass of 
the envelope increases. Ultimately, a maximum value of E is 
attained, beyond which no solutions exist# 

The equation describing small, radial adiabatic pulsations is 



6 . 




d 

dx^ 
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Tiq 
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^(1+b- J 


0, 


(6) 


Where 


(1-P) (7r]-8-2P) ^-3^ 


ri(8-7^ 


Ck 


n+iy' 


0,2 = 


c 


2tt n ^ 
Pery GM' 


(7) 


and § is the radial displacement (Stothers 1965) . vThe period- 
root-mean-density relation is simply 

Per =* (3 ttG/uj2)^* (8) 

At the surface V = so we require the following regularity 
condition, after normalizing the displacement, 

X = 1, 5=1, d?/dx = (uj^+4)/r;j^-2, (9) 

An eigenvalue iD^ is determined whenever the solution converges 
to 5 = 0 at the center (or inner boundary) . Rabinowitz (1957) 
has discussed methods of solution and the rapidity of convergence* 
Results of integrations of the structure and pulsation 
equations are presented in Table 1 over the permissible ranges of 
E and A. Pulsational eigenf requencies are given for the fundamental 
mode of pulsation only. The structures with no radiation pressure 
(A = 0) have already been treated by several authors (Osterbrock 
1953; Harm and Schwarzschild 1955). In a few cases their 
pulsational properties have also been calculated (Cowling 1934; 
Rabinowitz 1957) . It may be noted that the case E = 45.5 is 
simply the polytrope of index 3/2 with F = 5/3. The case 
E = 0 corresponds to an atmosphere of no mass, in which equations 
(4) are exact everywhere; however, the pulsation equation (6) may 
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still be solved in this case. 

A few structures including radiation pressure have been 
calculated by other authors. The structure with infinite mass 
(A = “) is simply the polytrope of index 3 with F = 4/3. The 
case E = 0 corresponds to an atmosphere of no mass with F = 4/3. 

A few non-degenerate structures have been calculated by Hayashi 
et al . (1962) . However , most of Table 1 is new. The maximum 
value of m anywhere occurs in the case A == 0 , E - 0 . 

The effect of radiation pressure is to deepen the convective 
envelope and to reduce the pulsational eigenf requency . Therefore, 
for the same envelope depth in mass fraction, more massive stars 
tend to have pulsational eigenf requencies quite comparable with 
those of less massive stars. Massive stars, however, can possess 
very deep convective envelopes during advanced stages of evolution 
(cf . references in Stothers 1968) . The physically realizable 
portion of Table 1 for massive stars lies between A - 10 and 
juUU . Tnese limTtl cbrrespoid, respectively, to stars of 
IOMq with an age-zero chemical composition and 60 Mq with a heavily 
me tal— enriched chemical compos ition . The possible values of Q are 
therefore confined to the range 

0.05 < Q < 0.10 day. 

It is somewhat surprising that the derived upper limit for Q is 
not considerably larger, in view of the large radiation pressure 
at high masses. 
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III. VARIABLE RED SUPERGIANTS 

The known variable M supergianta in galactic associations 
ore listed in Table 2. The spectral types and magnitudes are 
those adopted in the accompanying paper (Stothers 1968, Table 1) . 
The magnitudes selected as applicable are median magnitudes. The 
effective temperatures and bolometric corrections (for the various 
spectral subdivisions) have been taken from Johnson (1966) The 
radii can then be derived in the usual way from the effective 
temperatures and bolometric magnitudes. The periods are the ones 
listed in the General Catalogue of Variable Stars (Kukarkin, 
Parenago, Efremov, and Kliolopov 1958). The mass of an individual 
red supergiant has been assumed equal to the average (or, in some 
cases, the most frequent) mass of blue supergiants in the same 
association (Stothers 1968) . For each variable we have formed an 
observational Q value. 

According to the pulsation theory of Section II, these 
variables should have Q values in the'. range 0.05 Q ^ 0.07 day 
if they are undergoing radial pulsation in the fundamental mode. 
This range corresponds to the portion of Table 1 covered by the 
relevant A values for 14 to 25 M©. Since the convective envelopes 
in published models of red supergiants (cf. references in Stothers 
1968) are of intermediate depth, we can safely assume a uniform 
Q value. 
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Q =0.06 day. 

Five of the ten variables with listed periods in Table 2 
/ 

have Q values in reasonable accord with the pulsation hypothesis. 
Only one of the ten, RS Per, has a Q value which is significantly 
smaller ((corresponding approximately to a fifth overtone!). The 
other four have incompatibly large Q values. If it is assumed 
that the light variation in all cases is due to pulsation, _ then 
we must seek sources of possible error in our evaluation of Q for 
these stars. 


In terms of quantities directly used to determine Q, we have 

. -3/4 „ 3 




Per. 


The uncertainty associated with each of these quantities will be 


discussed in turn. 


Luminosity . The errors in the derived bolometric magnitudes 
are due to (1) observational errors and (2) association back-to- 
front differences and possible association non-membership. The 
observational error (Stothers 1968) is io.3 to ±0.7 mag., which 
corresponds to an error in Q of only 20 to 50 per cent. This is 
unimportant. In the outer I Per aggregate, T Per and S Per are 
very likely to be members of the background Perseus arm, lying at 
(m-M) ^ = 12.7 (Schild 1967). In the I Ori association, a Ori is 
probably a foreground member with My ~ -5.7 to -6.1 (Wilson 1959; 
Jenkins 1963) . For T Per, S Per, and a Ori, we have used the 
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adjusted luminosities to obtain the Q values listed in Table 2, 
rather than the luminosities in Stothers (1968) . The other stars 
questioned by Stothers — TV Gem, a Sco, and RS Per — are probably 
normally associated with their respective aggregates. Although 
the relatively late spectral type (M4.5) and therefore possibly 
low luminosity of RS Per suggest that it may be a, foreground 
star (larger Q value) , the star lies in the center of x Per 
against the sky and is almost surely a cluster member. 

Mass . In order to bring the" observed values of Q for a Sco, 
a Ori, T Per, and S Per into line with the pulsation hypothesis, 
the masses of these stars would have to be reduced to values 3 
to 100 times smaller than those listed in Table 2. Such a large 
amount of implied mass loss is inadmissible. Most of the stars 
would have radiative envelopes and appear blue in color, not red, 
apart from the unlikelihood that mass loss would affect certain 
stars of a given spectral type so much more than other stars of 
the same spectral type. Moreover, the, observed mass-loss rate 
for a Ori (Weymann 1962) is quite incapable of causing such a 
drastic amount of mass loss. In the anomalous case of RS Per, 
its present mass would have to be 25 times greater than its 
adopted mass to explain the discrepancy. This is also inadmissible 

Effective temperature. The scale of effective temperature 
as a function of spectral type is not known very well for M 
supergiants. However, by intercomparing the various determin- 



11 . 


ations of effective temperature for a Sco, a Ori , and Her 
(Johnson 1966 ) , We find that the error in the adopted value is 
probably less than i300°K. This causes a possible error in Q 
of only 30 per cent. In any case, the pulsationally discrepant 
stars are evenly distributed among the same spectral subdivisions 
as the concordant stars. 

Period . It is likely that virtually all the discrepancies 
in Q in Table 2 are due to errors in the listed periods. The 
long cycles and relatively small amplitudes (typically 1 mag. 
photographic) of variable M supergiants are very difficult to 
determine and are usually semiregular in nature. Multiperiodicity 
(or simply irregularity) is not unknown; for example, Palmer (1939) 
lists periods of 11, 250, and 2070 days for a Ori and 840, 1120, 
and 3360 days for S Per. 

With our adopted periods, no period-luminosity nor period- 
spectrum relation emerges from the data. Considering, however, 
only the five stars with Q values consistent with the pulsation 
hypothesis, we find both such relationships emerging (if YZ Per 
is assumed to be slightly in the foreground of the outer I Per 
aggregate) . Since the Q value is determined principally by 
the period and luminosity, the discovery of a period-luminosity 
relation for stars selected by Q value is not surprising. However, 
since our Q values are also dependent on spectral subdivision 
(through the effective temperature) , the discovery of a period- 
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spectrum relation could be further evidence in favor of the 
pulsation hypothesis. 

The reality of such a relation may be tested by adjusting 
the periods of all stars with "discrepant" Q values to fit an 
adopted Q = 0.06 day. The resulting "revised" periods are listed 
in the ninth column of Table 2. It is immediately clear that a 
period- spectrum relation is now supported by the data. This is 
basically a consequence of the luminosity-spectrum relation among 
the variables, the brighter being later in spectral type. It remains 
to be seen, of course, whether the observations permit such a 
readjustment of the periods to be made. In the case of a Ori, 
the revised period agrees well with one of the periods listed by 
Palmer (1939) . In the other cases (except for T Per) , the revised 
periods are nearly integral multiples or divisors of the original 
periods. Concerning the distances of a Ori, T Per, S Per, and 
YZ Per, our suspicions now seem to have been well founded since 
the revised distances conform well with the pulsation hypothesis. 
Finally, we mention the interferometric and radial- velocity 
studies of a Ori (Pease 1931; Sanford 1933) which suggest that 
this star undergoes a cyclical variation in radius. 

The relation between luminosity, spectrum, and period for the 
semiregular red variables is shown in Figure 1. It should be 
remarked that the visual absolute magnitudes of these stars are 
very similar (because of the effect of the bolometric correction) , 
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averaging <My> = -5.4. Since more than half the stars listed 
^in Stothers (1968* Table 1) with spectral type Ml and later are 
known variables, it is possible, and even likely, that the 
remaining stars are also variables of small amplitude. 

We conclude that the hypothesis of radial pulsation (in the 
fundamental mode) in the variable M supergiants tends to be 
supported by the available observational data. 
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Theoretical Convective Envelopes Including Radiation Pressure 
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FIGURE C;\PTION 


Fig. 1* “ Bolometric H-R diagram for the known variable M super- 
giants in galactic associations. The variables are 
labeled with their period in days, if known. Periods 
“revised" in accordance with the text are given in 
parentheses. 
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